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ABSTRACT 

Galaxy-scale strong gravitational lens systems are useful for a variety of astrophysical applications. 
However, their use is limited by the relatively small samples of lenses known to date. It is thus 
important to develop efficient ways to discover new systems both in present and forthcoming datasets. 
For future large high-resolution imaging surveys we anticipate an ever-growing need for efficiency 
and for independence from spectroscopic data. In this paper, we exploit the clustering of massive 
galaxies to perform a high efficiency imaging search for gravitational lenses. Our dataset comprises 
44 fields imaged by the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS), each of 
which is centered on a lens discovered by the Strong Lens ACS Survey (SLAGS). We compare four 
different search methods: 1) automated detection with the HST Archive Galaxy-scale Gravitational 
Lens Survey (HAGGLeS) robot, 2) examining cutout images of bright galaxies (BGs) after subtraction 
of a smooth galaxy light distribution, 3) examining the unsubtracted BG cutouts, and 4) performing 
a full- frame visual inspection of the ACS images. We compute purity and completeness and consider 
investigator time for the four algorithms, using the main SLAGS lenses as a testbed. The first and 
second algorithms perform the best. We present the four new lens systems discovered during this 
comprehensive search, as well as one other likely candidate. For each new lens we use the fundamental 
plane to estimate the lens velocity dispersion and predict, from the resulting lens geometry, the 
redshifts of the lensed sources. Two of these new systems are found in galaxy clusters, which include 
the SLAGS lenses in the two respective fields. Overall we find that the enhanced lens abundance 
(30lg^ lenses/degree^) is higher than expected for random fields (12^2 lenses/degree^ for the COSMOS 
survey). Additionally, we find that the gravitational lenses we detect are qualitatively different from 
those in the parent SLAGS sample: this imaging survey is largely probing higher-redshift, and lower- 
mass, early- type galaxies. 

Subject headings: gravitational lensing - techniques: miscellaneous - surveys 



1. INTRODUCTION 

Strong gravitational lensing - when the potential of a 
massive foreground object causes the formation of multi- 
ple images of a background source - is a powerful tool for 
cosmological and astrophysical research. Applications in- 
clude measuring the mass distributions of dark and lumi- 
nous matter, and measuring cosmologic al parameters via 
the lens geometry or abundance (see e.g. lKochane3l2QQ6l . 
for a review). A homogenous, well-understood sample of 
lenses is required for a statistically signific ant study, ne- 
cessitating large-scale, systematic surveys (jTurner et al.l 
[1983 lBolton"etlII [20061 : llnada et al.l[2008h . The - 200 
galaxy-scale lenses known today have been found through 
the numerous selection algorithms described in the next 
paragraph, and include a significant number of serendip- 
itous discoveries. Although this has enabled substantial 
progress, the number is still the limiting factor for many 
applications. 

Most searches so far have focused on either the source 
or lens population, employing a range of different strate- 
gies. Searches targeting potential sources have included 
looking for multiply-imaged radio so urces - as in the 
Cosmic L ens All-Sky Survey fCLASS: [MVers et al.l[2003l : 
Browne et al. 200 31) ~ a-nd exarnining known quasars (e.g. 
Maoz et al.l 119971 : iPindor et aLll2003[ ). The Sloan Digi- 
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tal Sky Survey (SPSS) Quasar Lens Search (|Oguri et al.l 
120061 : llnada et al.ll2008D selects lensed quasar candidates 
using two algorithms, one based on morphology and 
color for small-separation images, and one based only 
on color. Of the surveys targeting potential lenses, 
most have involved insp ectio n of high resolution ima ges. 
iRatnatunga et all (|l999f ) and lMoustakas et all (|2007[ ) se- 
lected candidates by eye from HST color images of the 
Extended Groth Strip (EGS) fields. Other authors have 
attempted to pre-select massive galaxies by their optical 
magnitude and color, and then examine the residuals af- 
ter galaxy subtraction. This was done in the Great Ob- 
servatories Origins De ep Survey (GOODS) lens search 
(|Fassnacht et al.l [20041 ) and t he Cosmic Evolut ion Sur- 
vey (CO SMOS) l ens se arch (|Faure et all 120081 ). Most 
recently, iJacksonI (|2008l ) extended this method to lower 
mass galaxies, but focused on viewing large- format ar- 
rays of unsubtracted col or galaxy cutout images. The 
Sloan Lens AGS Survey (|Bolton et al.ll2006l . [2008a[ ) can- 
didates were selected based on spectroscopic data indi- 
cating multiple redshifts in the spectrum of early-type 
(and hence massive) galaxies, then classified using this 
data and visual examination of high resolution HST im- 
ages (before and after lens galaxy subtraction). 

Recently, several algorithms for automated lens 
detection h ave b ee n developed. Tho se dev eloped 
by lAl ard (2006), "Seidel BartelmannI (|200/t ) and 
iKubo Deir An tonio (200j) look for arcs, a common 
feature in both weak and strong lensing on group and 
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cluste r scales. lEstrada et al.l (|2QQ7[ ) and lBelokurov et al] 
(|2QQ7l ) suggest mining databases for blue objects near 
potential lenses, since the most common sources are 
faint blue backgrou nd galaxies. The RingFinder 
(|Cabanac et al.| [2QQ?) applies the same logic to smaller 
image separation lenses, subtracting a rescaled red im- 
age from a blue cutout image to dig within the lens light 
distribution; it then analyzes the shapes and positions 
of the remaining residuals. Most re cently the HAGGLeS 
automated lens detection "robot" ^M arshall et al.ll2QQ8l ) 
attempts to model every object (typically selected to be 
Bright Red Galaxies or BRGs) as a gravitational lens, 
i.e. as a combination of background light from the source 
that is consistent with having been multiply-imaged, and 
residual foreground light from the lens galaxy. The re- 
sult is the robot's quantitative prediction of how a human 
would have classified the candidate. 

Fortunately, we are about to enter an era when orders 
of magnitude increases in the number of known lenses will 
be possible. In the near future, wide-field surveys such 
as the proposed Joint Dark Energy Mission (JDEM) and 
Euclid space missions would provide ~ 1 0^ — 10^ square 
degrees of high resolution imaging data ([Marshall et al.l 
I2OO8) . Automation will be needed to examine such large 
areas over manageable timescales, and to draw attention 
to those systems which have a higher probability of being 
lenses. We also must be prepared to proceed without the 
help of spectroscopic data. 

We describe here just such a survey: based purely on 
imaging data, and with a sufficient degree of automation. 
We compare the accuracy of four methods of searching 
for strong gravitational lenses. Presented in order of de- 
gree o f automation, they ar e: (1) using the HAGGLeS 
robot (iMarshall et al.ll2QQ8p. (2) examining; subtraction 
residuals (e.g. Faure et al. "2008^, (3) looking at galaxy 
cutouts (e.g. Fassnacht et al. 2004; Jackson 2008), and 
(4) performing a visu al inspection of entire fields (e.g. 
iMoustakas et al.ll2007f ). 

We additionally aim to make use of galaxy clustering 
in order to impr ove the efficienc y of our search, as sug- 
gested by .Fassnacht et al.l (|2006[ ) . The most likely lens- 
i ng galaxies are massive ellipticals with .3 < z < 1.3 
([Turner et al.lll984l : iFassnacht end] 120041 ): we may ex- 
pect that focusing on bright red galaxies (BRGs) will in- 
creas e the lens detection efficienc y (e.g. [Fassnacht et al.l 
120041 : iFaure et all l2008l : MarshalteLaD I2008D . The 
sources, meanwhile, are predominantly expected to be 
faint blue galax ies (FBGs) at redshifts at or above 1 
([Marshall et al.l l2Q05). Both of these types of objects 
(potential lenses and potential sources) are clustered, 
BRGs more strongly that FBGs. We therefore ex- 
pect that strong gravitational lenses are also clustered; 
we anticipate finding more new lenses by looking near 
known ones than we would otherwise. For example, 
IFassnacht et al.l ([20061 ) presents just such an occurrence: 
the researchers discovered two additional lens candidates 
less than 40 arc-sec from the known lens B1608+656. 

Therefore, as our dataset we use a subset of the SLAGS 
HST/ACS fields. The SLAGS surv ey has discovered 70 
definite galaxy-galaxy strong lenses ([Bolton et al.ll20Q8al ) 
to date. These definite lenses have clearly identifiable 
lenses arcs or multiple images in addition to spectro- 
scopic data; due to the rigorous requirements for this 
classification, we consider all to be confirmed lenses. Of 



these, 63 are well- modeled by a singular isothermal el- 
lipsoid (SIE) and have lens and source redshifts along 
with F814W photometry for the lens gala xy; most also 
have m easured stellar velocity dispersions ([Bolton et al.l 
'2008a'). This makes the SLAGS lenses the largest ho- 
mogenous sample of strong lenses to date. Each of the 
fields we use is centered on one of the known SLAGS 
lenses (hereafter the "main lens"); we thus make use of 
gravitational lens clustering in every field. The majority 
of our fields were observed in just one filter, meaning that 
we by necessity pre-select bright galaxies (BGs) rather 
than BRGs. For the automated portion of our searches, 
we use the HAGGLeS robot which, in addition to per- 
forming a quantitative lens classification, creates a useful 
database of galaxy cutouts and residuals. 

The organization of this paper is as follows. In Sec- 
tion [2] we present the 44 AGS fields in our dataset, and in 
Section [3] we outline the four search procedures we used 
to find gravitational lenses in these images. To assess 
the accuracy of these methods, we discuss each proce- 
dure's performance on the main lenses in Section [H We 
then present our four new definite strong lensing systems 
and one likely candidate, and investigate their physical 
properties and environments, in Section O After some 
discussion of our methods and results in Section [H we 
conclude in Section [71 All magnitudes are in the AB sys- 
tem; we assume a fiat cosmology with 1^^=0.3, 1]a=0.7, 
and Ho=100/i kms~^Mpc~^, where h=0.7 when neces- 
sary. 

2. SAMPLE SELECTION AND OBJECT DETECTION 
We use a subsample of the ACS F814W fields in- 
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sample selection and data reduction and analysis of the 
SLAGS data is given in SLAGS-I, -IV and -V and wih 
not be repeated here. We select the subset of 44 definite 
lenses whose ACS observations meet the following two 
requirements: exposure time >1000 seconds through 
F814W, and lens classification "definite." Due to the 
demise of ACS, two filters are not available across all 
fields; of the AGS filters used, F814W is the most widely 
available in the SLACS fields. 14 fields were additionally 
imaged through F555W. We require long exposure time 
to guarantee the ability to produce deep, high quality, 
cosmic ray- free reduced images. 

Our final sample comprises 44 uniformly-observed ACS 
fields, each centered on a SLAGS gravitational lens 
(a "main lens"). Each field covers approximately 11 
arcmin^; thus our total survey area is 0.134 square de- 
grees. 

Most lens-finding methods involve examining bright 
galaxies (BGs ), as being the most li kely to be lensing 
galaxies (e.g. 
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Fassnacht et"all 120041 : iFaure et all l2008l : 
I2OO8I ). Here this selection is done with 



magnitude and size cuts on catalogs made with the SEx- 
tr actor ^ software tool; only a third of our fields have 
multi-filter ACS data so in order to maintain consistency 
we do not use color selection even when it is possible. We 
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choose an apparent F814W magnitude (mF8i4w; using 
Kron-like elliptical apertures) limit of 22. Other limits 
were tried, but 22 was found empirically to be the best 
balance between efficiency and the selection of all bright 
galaxies: with mF8i4W < 22, a total of 2399 BGs are 
selected, but with mF8i4W < 23, that number jumps to 
5439 without adding significantly to the actual number of 
morphologically early- type galaxies selected. The mini- 
mum size selected is FWHM=0'.'215, in order to reject 
most stars. We additionally use the SExtractor "flags" 
parameter to limit the remaining contaminants. 

3. SEARCH METHODS 

We use four distinct methods to search for gravita- 
tional lenses. They are presented in order of degree of 
automation, where the first has the most pre-selection 
and analysis, and the last none: 

• Procedure 1: Using output from the HAGGLeS 
robot 

• Procedure 2: Inspection of subtraction residuals of 
bright galaxies 

• Procedure 3: Inspection of cutout images of bright 
galaxies 

• Procedure 4: Full-frame visual inspection of the 
ACS images 

In general we expect our searches to be sensitive to 
lensing events similar to, but fainter than, the SLAGS 
main lenses. With the BG selection in Procedures 1-3, 
and the innate limits of the human eye in Procedure 4, 
we are biased towards lensing events around large, lumi- 
nous galaxies. However, empirically the main lenses are 
the brightest objects (other than saturated stars) in the 
fields, so we expect any other lenses we find to be, by ne- 
cessity, fainter. The main lenses also favor high magnifi- 
cation configurations (Einstein rings) due presumably to 
their spectroscopic selection; we expect to be somewhat 
preferential in this regard as well, on the grounds that 
these distinctive cases will be easiest to identify visually. 

Before continuing, we first consider the question "what 
is the probability that this object is a gravitational lens?" 
The classification of lens candidates has varied among 
the major surveys targeting potential lenses, as has the 
data on which th ese classifications are based. The GOS- 
MOS lens search (|Faure et al.ll2QQ8l ) sample is subdivided 
into "best systems" , which are deemed to have a greater 
probability of being lenses than the remaining objects. 
GOODS (Fassnacht et al. 200^ s elected candida tes and 
voted to choose "top candidates." iJacksonI (|2QQ8l ) groups 
potential lenses into three categories: very likely or cer- 
tain lenses, possibl e or probable lenses , and not-lenses. 
The SLAGS survey (iBolt on et "aIll2QQ8af ) also groups can- 
didates into three categories (definite lenses, probable 
lenses, and inconclusive/not lenses). We note that the 
classification of definite or probable lenses is not the same 
across surveys, and therefore one has to be careful to 
impose similar quality criteria when comparing inferred 
density of lens galaxies. 

In this work we follow [Marshall et al.l (|2QQ 8h and em- 
ploy a 4-point subjective classification scheme, outlined 
in Table [Tj the classification parameter H may range 



TABLE 1 
Lens Classification Systems 



Human class H 


Robot class Hr 


Definition 


3 


>2.5 


Definitely a lens 


2 


1.5 to 2.5 


Probably a lens 


1 


0.5 to 1.5 


Possibly a lens 





<0.5 


Definitely not a lens 



from 3 (definite lenses), through 2 (probable lenses) and 
1 (possible lenses), to = (definitely not lenses). Each 
of the 44 fields in our sample contains a confirmed lens 
from the SLACS survey, where here the grade is based 
on all available data, notably the galaxy spectrum (in- 
cluding anomalous high redshift emission lines) and clear 
lensing-consistent residuals after lens galaxy light sub- 
traction in all procured filters. We thus assign the clas- 
sification parameter = 3 for each main lens. 

Humans (after some training) are adept at identifying 
gravitational lenses by eye, using high resolution imag- 
ing data alone; we effectively make an internal model for 
the lens and optimize it. However, the amount of infor- 
mation available during each lens search with which to 
do this may vary: in any given search we always lack 
additional understanding that would aid us in identify- 
ing lenses. When carrying out the four different search 
procedures, we therefore assign each system a value Hi 
(where i = 1, 2, 3 or 4): each of these is our best guess as 
to the classification that a trained human would have 
assigned the system, if they had been given only the 
data presented in the i^^ procedure. (This results in four 
"re-classifications" for each main lens.) We consider the 
"true" classification H to be the value we give an object 
when taking into account all data. 

3.1. Procedure 1: HAGGLeS robot 

In Procedure 1, we use the HAGGLeS robot, 
the automated lens detection program developed by 
[Marshall et al.l (|2QQ8l ) , to identify samples of lens candi- 
dates prior to visual inspection. The robot treats every 
object as if it were a lens, models it, and then calculates 
how well the lens hypothesis works. It extracts and uses 
6 arc-second square cutout images of each object, and 
then focuses on the residuals made by subtracting off an 
elliptically symmetric Moffat profile model for the puta- 
tive lens galaxy light. The gravitational potential of the 
lensing galaxy is assumed to be sufficiently well-described 
by a singular isothermal sphere (SIS) plus external shear; 
this model is fitted to the residual image. Disk-like fea- 
tures and bright neighboring objects are masked before 
fitting. When multi-filter data is available it is used, and 
although the robot does not rely on color, with it per - 
formance is expected to improve (jMarshah et al.ll2QQ8[ ). 
Based on the results of the modeling process, the robot 
calculates and assigns a value of Hr to each object, where 
Hr is the robot's estimate of the classification H a human 
would have given the system. We summa rize the human 
and r obot classifications in Table [T] (after [Marshall et al.l 

Eool). 

The robot may have one of two characters, reflecting 
the prior probability that an object is a lens. The "re- 
alistic" character robot might expect 0.1% of objects to 
be lenses, and the "optimistic" character robot might 
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TABLE 2 

Prior probability distributions for robot characters 



Character 


Ft{H=0) 


Ft{H=1) 


Ft{H=2) 


Ft{H=3) 


Realistic 


0.900 


0.080 


0.019 


0.001 


Optimistic 


0.050 


0.100 


0.250 


0.600 



expect 60% to be (|Marshall et al.ll2QQ8D . The realistic 
robot is approximately in accordance with current es- 
timates of th e fractio n of strongly lensed galaxies (e.^. 
i Bolton et aT] [2QQ6: M oustakas et alJl2QQ7l : lMarshan et all 
|2008'), while the optimistic robot has the advantage of be- 
ing most inclusive by giving higher values of Hr. These 
two p rior probability distr ib utions are given in Ta ble [2] 
f after iMarshah et alJ[2QQ8h . [Marshall et all (|2QQ8r ) note 
that while the realistic robot produces lens samples of 
high purity - and correspondingly little need for hu- 
man inspection - a highly complete search requires, at 
present, a more optimistic robot. 

We ran both the optimistic and realistic robots in order 
to verify this last claim, but chose to use the optimistic 
robot to search for new lenses, as our sample is small 
and we are most interested in finding all lenses present 
in our sample. When comparing the performances of 
both robots, we refer to the optimistic and realistic robot 
classifications as H^^^ and H^^^^ , respectively. The main 
goal of the comparison is to gain insight into the pop- 
ulations of objects selected, and thus help improve the 
automated part of this procedure for future surveys. 

The output from the robot includes an individualized 
webpage for each object, containing the cutout, subtrac- 
tion residual, masked residual, and estimation of Ein- 
stein radius. This page also displays the minimal source 
able to produce the observed configuration, and the pre- 
dicted image-pla ne residual reconstruc ted from this min- 
imal source fsee iMarshah et al.l [20081 for more details). 
We viewed the pages of objects for which H^^^ > 1.5 
(indicating optimistically probable and definite lenses), 
and gave each object a human classification i^i, using 
the same data that was available to the robot (namely, 
the cutout image and subtraction residual image), and 
the lensmodel outputs produced by the robot. 

3.2. Procedure 2: Examining BG residuals 

Procedure 2 makes use only of the subtraction residu- 
als produced by the robot. We examine the residuals of 
every BG initially selected in the fields; all of the by 
6" subtraction residual cutouts for a field are displayed 
in a grid for rapid viewing. When multi-filter data is 
available, the residuals are shown in color. Looking at 
the 44 grids, we note all objects of interest; to each of 
these objects we assign a human classification i^2, based 
solely on the subtraction residual data. 

3.3. Procedure 3: Examining BG cutouts 

In this method, we inspect each of the by cutouts 
of the BGs provided by the HAGGLeS robot. We again 
display the cutout images in a grid for rapid viewing, 
making t his pr ocedure very similar to that employed by 
IJacksonl (|2008t ). However, the stretch of the cutouts' im- 
age display is fixed for al l objects at a level app ropriate 
to the EGS-type lenses (jMoustakas et al.ll2007D . This 



makes this method particularly sensitive to faint lens- 
ing events, and somewhat insensitive to lensing by very 
bright galaxies. We assign objects of interest a human 
classification, this time i^s, using only the galaxy cutout 
on display. 

3.4. Procedure 4-' Full- frame visual inspection 

During the full-frame visual inspection, each of the 44 
fields' F814W ACS images is viewed with ds9.^ We ini- 
tially set image parameters in each field so as to be sen- 
sitive to the lensing events we expect to find: similar to, 
but fainter than, the main lenses. We do this by setting 
the intensity stretch such that the main lens is slightly 
over-saturated; from this display we select objects of in- 
terest. For each object, we adjust the scale limits and 
intensity stretch such that the potential lensing features 
are most apparent, then decide on a classification. Pro- 
cedure 4 differs from the previous three most notably in 
that all galaxies in the field - not just BGs - are exam- 
ined, and that their viewing parameters are set individu- 
ally. The human classification parameter assigned in this 
procedure is denoted H^. 

3.5. Summary of Results 

During each procedure, any potential gravitational lens 
was marked and assigned a human classification based 
only on the information available during the search pro- 
cedure by which it was found. Following the completion 
of all searches, all available information was considered 
holistically and objects were assigned a final human clas- 
sification H. The classification is based primarily on our 
ability to recognize a typical lens geometry in the ob- 
served morphology (in both the cutout and subtraction 
residual), on colors and surface brightnesses consistent 
with lensing, and on the robot's ability to model the ob- 
ject. In total we discovered four new objects with i7 = 3, 
and one object with H = 2. These new gravitational 
lenses are presented in Section [5l 

4. METHOD ACCURACY 

Before investigating the properties of our newly- 
discovered lenses, we first discuss the performance of the 
four search procedures in terms of both their ability to 
find SLACS-type lenses (low redshift, high magnifica- 
tion), and their contribution to the finding of the new 
gravitational lens systems. For each method, we calcu- 
late the purity and completeness of the samples selected, 
considering our five best systems and the SLAGS main 
lenses together as lens systems to be recovered. We de- 
fine purity as the percent of selected objects (i.e. hav- 
ing classification parameter Hi above some threshold), 
that actually have final human classification H greater 
than the same threshold. Only the SLAGS main lenses 
and our 5 best systems have H > 1.5, and all but one 
of these has H = 3. Similarly, completeness is defined 
as the percent of the objects with final human class H 
greater than some threshold that were given procedure 
classification Hi greater than the same threshold. For 
Procedure 1, any object that was not examined because 
H^P^ < 1.5 is considered to have Hi = 0. Purity and 
completeness were calculated for each procedure, and for 

^ http: / / hea-www.harvard.edu /RD /ds9 / 
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TABLE 3 

Purity and completeness for each method 



Procedure 


Hi cut 


Purity(7oj 


Completeness (To J 


1 


Hi > 2.5 


100 


69 




Hi > i.o 


iUU 


on 


2 


^ 9 
112 Zi.o 


100 


77 




H2 > 1.5 


100 


90 


3 


H3 > 2.5 


100 


4 




H3 > 1.5 


93 


22 


4 


H4 > 2.5 


100 


81 




H4 > 1.5 


100 


92 


Rea. robot 


if;^"^ > 2.5 


40 


8 




if^^"^ > 1.5 


8 


41 


Opt. robot 


> 2.5 


5 


27 




> 1.5 


3 


82 



both the realistic and optimistic robots. These statis- 
tics are summarized in Table [3l errors are of order a few 
percent. 

Note that during each search procedure, the main 
SLAGS lenses themselves ae classified, resulting in values 
of i^i, i^2, ^3, and H4 for each. A lens-by-lens compar- 
ison of classification values is available in the Appendix, 
while data on the total numbers of main lenses and new 
lens recovered by each search procedure are presented 
in Table [H In this table we also include the amount of 
time spent per field for each method, and include the 
performance of the HAGGLeS robot (with no human in- 
spection) for comparison. 

In the subsections that follow, we briefly discuss these 
results and their implications for optimizing lens discov- 
ery in high resolution imaging surveys. In the Appendix 
we show, for reference, some of the images that com- 
prise the input data to the first three search procedures, 
namely, cutout images of the main lenses, and of the 
i?;^^^ > 2.5 candidates. 

4.1. Procedure 1 

Inspecting all objects classified by the optimistic robot 
as rn^'P^ > 1.5, we classified three out of the five new sys- 
tems as iYi = 3, and one system as i^i =2. The final 
system had been assigned a robot classification parame- 
ter of H^P^ < 1.5 and so was not examined. This failure 
is most likely due to the non-lensed light remaining af- 
ter the Moffat subtraction (see the Appendix for more 
illustration of this). 

Indeed, for the main lenses, we found that the un- 
masked disk features are the most common cause of 
robot detection failures: we show the subtraction resid- 
uals and reconstructed images in the Appendix. The 
robot may fail to find a suitable model entirely in some 
extreme cases, such as when considering objects with 
very strong disk components or bright companion galax- 
ies (main lenses SDSSJ1029+0420, SDSJ1103-f5322 and 
SDSSJ1416+5136). We find overah that 8 of the main 
lenses, including the three cases just mentioned, are 
not detected at H^^^ > 1.5 by the robot. Interest- 
ingly, one of the 36 robot-detected systems, and three 
of the 8 robot-rejected systems were human-classified 
as i^i < 1.5: given just the output from the robot, 
these lenses (SDSSJ1213+6708, and SDSSJ1016+3859 
SDSSJ1029+0420, and SDSSJ1032+5322) would not 
have been identified as lenses by a human inspector. 

The lens galaxy subtraction also seems to be the major 



cause of false detections: ring-like and disk-like features 
left over from the Moffat profile subtraction and then 
incompletely masked can be wrongly interpreted by the 
robot as lenses. Examples of such false positives can be 
seen in the 10 objects classified by the realistic robot as 
ij;^^^ > 1.5 (in the Appendix). 

Of the 36 robot-detected main lenses, only 20 were 
classified by the optimistic robot as "definite lenses" 
{H^^^ > 2.5). The reasons for these mis-classifications 
are a little more subtle. The probability density func- 
tions (PDFs) used to calcul ate the Hr valu e s for this re- 
search were determined by [Marshall et al.l (|2008l ) based 
on a training set of EGS non-lenses and simulated lenses. 
By overlaying the robot model output for the SLAGS 
main lenses on these PDFs, we can gain insight into the 
cause of of the robot's mis-classifications. In Figure [T] 
we show Pr(d|i^ = 3), and overlay the d- values for the 
main lenses (larger data points). Here we clearly see that 
the PDF (contours approximating the density of smaller, 
training set points) is not optimized for SLAGS-type (low 
redshift, high magnification) lenses: to the robot, source 
magnitudes are surprisingly bright and arcs are of un- 
usual thickness. These differences are to be expected, but 
we note that the lenses we should expect to find with the 
HAGGLeS robot would therefore be more similar to the 
EGS lenses in terms of apparent magnitude and geomet- 
ric configuration. In a future wide-field search, where 
more SLAGS-type lenses may be present, it would be 
prudent to retrain the robot on a wider variety of lenses 
- and perhaps on the SLAGS sample itself. 

Having detected a SLAGS main lens, how accurately 
does the HAGGLeS robot model it? As an illustration of 
the robotic models' performance, we compare the robot- 
calculated Einstein radii (^e) to those determined in the 
SLAGS papers (Figure [2|). The 12 open data points are 
main lenses that were not well- modeled by the robot; af- 
ter removing these we find an rms scatter of 12%. For 
60% of the main lenses, then, the robot finds not only a 
successful lens model, but one that is in very good quan- 
titative agreement with that inferred during the SLAGS 
project. 

4.2. Procedure 2 

When examining the BG subtraction residual images, 
we classified all four of the new H = 3 systems as 
H2 > 1.5 or above; the H = 2 system was missed. Ein- 
stein rings are particularly easy to identify, resulting in 
a high proportion of the main lenses found. However, 
we occasionally lost some image context: some residu- 
als are more easily identified as being caused by lensing 
when viewed along side of the lens galaxy cutout. When 
looking at both cutouts and subtraction residuals, we are 
able to ask the question "Is this structure a part of the 
lensing galaxy or is it unique?" However, since sources 
are typically bluer than the lens galaxy, color residuals 
can aid us in answering this question: we note that all 
but one main lens was classified as = 3 when multiple 
filters were available. 

4.3. Procedure 3 

We found that Procedure 3 was not an appropriate 
method to use when looking for occurrences like the 
SLAGS lenses. This is predominantly an issue of im- 
age intensity scale and dynamic range. The scales used 
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TABLE 4 

Summary of results for each method 



Procedure 



Rea. robot Opt. robot 



Main lenses Hi > 1.5 
Main lenses Hi > 2.5 
New lenses Hi > 1.5 
New lenses Hi > 2.5 
Typical candidates 

Time per field (min) 
Time per degree^ (hours) 



20 
4 





36 
10 
4 
3 



35 
30 

4 

3 

arcs 

multiple images 
2-4 
10-20 



39 
34 

5 

3 

spirals, arcs 
multiple images 
1-2 
5-10 



2 
3 


blue images 
companions 
1-2 
5-10 



42 
36 

3 

3 

spirals, arcs 
companions 
20-40 
100-200 




Fig. 1. — Pr(d|/f=3) derived from the robot outputs of the EGS training set (smaller points), overlain with the outputs 
of the SLAGS main lenses (larger points). Points c orrespond to objects with Hi=3 and the contours in this and the 
subsequent two figures are 68% and 95% CL.. After [Marshall et aP (jboOSD . 



Finding lenses near lenses 



7 




0.5 1 1.5 2 



(SLACS) 

Fig. 2. — Comparison between the SLACS-measured Einstein 
radii (SLACS-V) and those calculated by the HAGGLeS robot. 
The systems for which robotic lens modeling failed are represented 
by filled points. The one-to-one relation is also shown to guide the 
eye; we find an rms scatter of 12%. Errors on SLACS Einstein 
Radii are 2% (SLACS-V). 

by the robot are set for fainter lens galaxies than those 
in the SLACS sample, and as such the main lenses often 
appear saturated: the lensed features were completely 
washed out in all but 8 cases. We consequently classified 
only two objects as = 3, and this is the only procedure 
in which we falsely identified an object as a lens. This 
procedure was found to be most effective when lensed 
images are blue against a red galaxy: five of the eight 
main lenses given Hs > 1.5 were in color fie l ds. T he is- 
sue of intensity scaling was noted bv lJacksonl (|2QQ8l ). who 
attempted to optimize the intensity scale on an object by 
object basis, and hence per formed ra ther better in terms 
of "true" lenses recove red (iJackson recov ered ~ 50% of 
the lenses identified bv lFaure et al.ll2QQ8l ). However, the 
large dynamic range in surface brightness of the lens and 
source galaxies makes this a very difficult task. 

4.4. Procedure 4 

In Procedure 4 we found three out of four new H = 3 
lenses, and only two of the main lenses were not classiffed 
as H4 > 1.5. The greatest strength of this method lies in 
the ability to vary the intensity scale of the image so as 
to make apparent any lensed structures. It is effectively 
the same as looking at both a cutout and a subtraction 
residual: lensing can be viewed in the context of the 
host galaxy. Its weaknesses, however, are in the time 
required and the quantity of data presented: it is difficult 
to conduct a thorough, mistake-free search when every 
object must be examined and the time required is ^30 
minutes per field. 

4.5. Considerations and comparisons 

We find that procedures one and two are nearly iden- 
tical in their ability to find new lenses, while the latter 
performed marginally better on the main lenses since all 
the objects were examined. Procedure 3 we find unsuited 
to the task at hand; looking at the robot output (Pro- 
cedure 1) or subtraction residuals (Procedure 2) largely 



circumvents the dynamic range problems of image dis- 
play, and produces superior results. Procedure 4 is the 
most inclusive of the four procedures as all galaxies (not 
just bright galaxies) are examined; it is also the only 
procedure in which we would have a chance of finding 
a "dark lens." However, we found two additional lens- 
ing systems when we pre-selected bright galaxies and our 
search took significantly less time. We hence find that 
looking only at the BG subsample does not decrease the 
completeness of a survey, and could in fact improve it 
due to a lower error-rate. 

Taking into consideration efficiency, completeness, and 
purity, we recommend the use of Procedure 1 and Proce- 
dure 2. There will be cases, such as with eight of the 
main lenses and one of our new lenses, where a lens 
will be missed by the current robot and thus by the 
human following Procedure 1. There are also certain 
lenses, particularly very unusual ones, that the current 
robot will miss due to its inability to model it. For ex- 
ample, the naked cusp configuration that often occurs 
with edge-on spirals produces three images blended into 
an arc but no counter-image; without a counter image 
the current robot will classify the system as a class 0. 
Additionally, when the environment plays a strong role 
- in over-dense environments for instance - the simple 
SIS+external shear model used now may not be suffi- 
cient. It therefore may be preferential to use a combina- 
tion of procedures 1 and 2, inspecting all objects modeled 
by the robot regardless of the modeling outcome. 

5. RESULTS 

We have discovered four new definite gravitational 
lenses and one promising lens candidate. These five sys- 
tems were each found in at least one of the four separate 
search procedures. Having taken into consideration all 
available data-including spectroscopy for one system-we 
assigned four objects true human classification H = 3 
(our four new lenses), and only one H = 2 (our best can- 
didate). The five systems were also classified during each 
of the procedures. Classifications for potential lenses are 
referred to as i^i, i^2, ^3, or H4 according to procedure 
number; these values may differ from H. We present the 
systems along with their classifications in Table [5l For 
simplicity, each lens is given a short name that will be 
used in the remainder of this paper. 

5.1. Improved Robotic Lens Models 

As shown in Section [4j and noted by [Marshall et al.l 
(|2QQ8l ). the lens model parameters returned by the HAG- 
GLeS robot are not always accurate. We find that the 
principal cause of robot modeling (and indeed classifica- 
tion) error is insufficient lens light subtraction. Disks and 
irregular profile slopes both give rise to significant sym- 
metrical residuals that confuse the robot. For the small 
number of high quality lens candidates identified in sec- 
tion [5l we can solve this problem on a case by case basis, 
and provide the robot with cleaner images to model, and 
thus produce more accurate estimates of the lens candi- 
dates' Einstein radii. 

We improved the lens galaxy light subtraction with 
the fiexible B-spli ne fitting approach developed by 
iBolton et al.l ([2006) , after first masking out the objects 
in the image identified as candidate lensed arcs. This 
procedure leaves sharper lensed image residuals. We 
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TABLE 5 

Classification of new candidates by method 



Short name 


Object name 


H 


Hi 


H2 


H3 


H4 


Tjreal 
11 J. 


hT 


Danny 


HST J073736.40+321540.3 


3 


3 


3 


2 


3 


0.9 


1.9 


Frenchie 


HST J114331.46-014508.0 


3 




2 


2 





0.1 


1.1 


Kenickie 


HST J121346.57+670833.3 


3 


3 


3 


2 


3 


0.7 


2.9 


Sandy 


HST J143001. 28+410440. 8 


3 


3 


3 


1 


3 


0.2 


2.9 


Rizzo 


HST J110307.14+532042.6 


2 


2 


2 


1 





0.5 


2.7 



Frenchie was not examined during procedure 1 because < 1.5 for this system. 



then set all undetec ted pixels to zero as described in 
[Marshall et alJ (|2QQ8[ ), but also at this stage masked out 
all the features in the B-spline residual map not iden- 
tified as lensed arcs. We then re-ran the lens-modeling 
part of the HAGGLeS robot, and took as our final es- 
timated Einstein radius the position of the peak of the 
G aussian fit to the sour ce plane fiux curve, as described 
in [Marshall et alJ (|2QQ8l ). The resulting lens models and 
their parameters are shown in Figure [3] and Table [6l 

We can use our results from Section [4m to estimate the 
accuracy of the Einstein Radii measured by the robot. 
As shown in Figure [21 we found that - for well modeled 
lenses - robot- and SLACS-measured Einstein radii agree 
to within 10%. However, this can be considered an up- 
per limit to our true uncertainty, since we use improved 
subtractions for the new systems. In practice, the robot's 
estimates using improved subtractions are as good as the 
relatively simple models allow them to be. To be conser- 
vative we adopt an error of 5% on the robot's Einstein 
Radii. 

5.2. Photometry 

Available from SDSS^ are photometric redshifts and 
apparent lens galaxy magnitudes; with the exception of 
one (Frenchie), no spectra are available. To supplement 
the S PSS data, we use the galf it software (jPeng et al. 
[2QQ2f ) to fit de Vaucouleurs models to the Hubble F814W 
data and derive apparent magnitude and circularized 
effective radii, listed in Table [6l After correcting for 
Galactic extinction, we calculate the rest frame V-band 
magnitude from the observed F814W, a conversion for 
which there is very little scatter between different spec- 
tral types. Details on surface photometry a nd K- color 
correc tions can be found in the paper by [Treu et al] 
([2OOII). For the main SLAGS lenses we use data from 
[Bolton et aP (2008a). 

5.3. Velocity dispersion and source redshift 

In this Section we use available spectroscopy and pho- 
tometry to estimate the velocity dispersion and source 
redshifts of the new lens systems, using two physically 
motivated assumptions: i) early- typ e lens galaxies lie on 
the fundamental plane (hereafter FP[Dre ssler et all[l987 : 
DiorgQvski Davis[fT987[ : [Treu et al.[ [2b06: Bo lton et all 
2008b[ ): ii) the ratio between stellar velocity dispersion 
and that of the best fitting SIS is approximately con- 
stant. 

We use the photometric redshift and evolution- 
corrected V band luminosities to estimate the central 
stellar velocity dispersion predicted by the FP relation: 

logi^e = a log CTe, 2 + 61og/e + C, (1) 
^ http://cas.sdss.org/dr6/en/ 



where (Je,2 is the stellar velocity dispersion corrected to 
an aperture of radius half the effective radius (in units 
of 100 kms~^), le is the effective surface brightness in 
units of 10^ L0 kpc~^, and Re is the effective radius 
in kpc. We adopt the coefficients a = 1.28 b = —0.77 
c = -0.09 derived in SLACS-V for the SLAGS sample. 
The intrinsic scatter of the fundamental plane dominate 
the uncertainty on the estimated cre,2 (0.05 dex). Gen- 
tral stellar velocity dispersions are obtained from <Je,2 
using the standard correction described in SLAGS-V. 

Additionally, SLAGS-IV and -V found that a* is corre- 
lated with (JsiE, the velocity dispersion that best fits the 
model of the lens as a singular isothermal ellipse (SIE). 
For the SLAGS sample: 



(^*/^tsie) = 1.02 ±0.01. 



(2) 



We may also calculate asm, assuming the Einstein ra- 
dius 0E, lens redshift, and source redshift are known: 



CTSIE 



47r Dds ' 



(3) 



where Dg and Dds are, respectively, angular-diameter 
distances to the source galaxy and between the lens and 
source galaxies, and 6^ is given in radians. 

We combine Equations [1] to [3] to determine our best 
estimate of the source redshift. To obtain the posterior 
probability distribution function of Zs , we assume that a* 
is log-normally distributed with scatter 0.06 dex, which 
is dominated by the intrinsic scatter of the FP and of 
Eq[2l We adopt priors appropriate for the source pop- 
ulation. For the newly discovered - imaging selected - 
lenses, we adopt as prior the redshift distribution of faint 
galaxies in single orbit AGS-F814W data as measured 
by the GOSMOS survey. For the main SLAGS lenses 
we adopt the same prior, but truncated at Zg < 1.5, i.e. 
the highest redshift where [O II] is still visible within the 
observed wavelength range covered by the SDSS spectro- 
graph used for discovery. The results change very little 
if a uniform prior is adopted instead. 

The result of this calculation for the main lenses can 
be seen in Figures [H and [5l The former compares the 
estimated source redshifts to the known source redshifts 
for the main SLAGS sample. This sanity check indicates 
that our procedure is unbiased and that the scatter is 
consistent with the estimated error bars. The latter fig- 
ure shows the posterior probability distribution function 
for the source redshifts of the newly discovered lenses. 
As expected, the posterior is asymmetric with a tail to 
high-z due to to the strong dependency on the ratio of 
angular diameter distances on the source redshift, when 
it approaches the lens redshift. 

The estimated stellar velocity dispersions and source 
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TABLE 6 

Relevant lensing data for new lenses 



Short name 


Separation 




F814W 


Re 




Zs 




<^*,FP 


Morphology 




(arcsec) 


(arcsec) 




(arcsec) 






(km s-1) 


(km s-i) 




Danny- 


107.9±0.2 


0.65±0.03 


18.75 


0.61 


0.34±0.03 


0.48+0 02 




296±41 


SO 


Main lens 




1.00 


17.04 


2.82 


0.3223 


0.5812 


358±17 


313±43 


E/SO 


Frenchie 


46.8±0.2 


1.45±0.07 


16.10 


1.79 


0.104 


'^•^-'^-0.56 


248±17 


226±31 


E 


Main lens 




1.02 


14.96 


4.80 


0.106 


0.4019 


279±13 


269±37 


E/SO 


Kenickie 


35.2±0.2 


0.65±0.03 


19.99 


0.43 


0.33±0.08 


1 44+^-^^ 

-•-•^^-0.57 




162±22 


SO/Sa 


Main lens 




1.42 


15.60 


3.23 


0.123 


0.6402 


308±15 


253±35 


E/SO 


Sandy 


82.7±0.2 


1.31±0.07 


18.27 


1.55 


0.32±0.01 


l-^«_0.56 




233±32 


E 


Main lens 




1.52 


16.87 


2.55 


0.285 


0.5753 


343±32 


305±42 


Sa/S 


Rizzo 


106.0±0.2 


1.24±0.06 


19.47 


0.46 


0.41±0.03 


1 -19+0.12 
-^•-^^-0.37 




279±39 


E 


Main lens 




1.68 


16.43 


1.95 


0.158 


0.7353 


211±12 


247±34 


undetermined 



Notes: F814W magnitudes are not corrected for Galactic Extinction. Uncertainties are typically 0.03 on 
F814W total magnitudes, 5% on Re, and 2% on SLAGS Einstein Radii (iBolton et al.ll2008ar) . 



redshifts are given in Table [6l We include the spectro- 
scopically measured cr^ when available, for comparison. 

6. DISCUSSION 

We begin this section by summarizing the arguments 
that led to our classification of the five newly discovered 
lens systems in Section 16.11 We then discuss the envi- 
ronments in which the lenses are found in Section \6.2\ 
and in Section [631 we compare the properties of the new 
systems to the SLAGS main lenses. 

6.1. Validity and Properties of Candidates 

Our final classification is entirely imaging-based (ex- 
cept in Frenchie's case, where an SDSS spectrum was 
available). We thus rely entirely on lens geometry, the 
appearance of the subtraction residuals, and the robot's 
ability to model the lens light. For this reason our stan- 
dards are high: we require clearly-identified multiple im- 
ages in all cases, and a straightforwardly- modeled im- 
age configuration. When available, colors are used to 
strengthen the case. As an additional sanity check, we 
note that the Einstein Radii are consistent with those ex- 
pected from the FP and a simple SIE model, for sensible 
values of source redshifts (see Table [6]). We now discuss 
each case individually. 

Danny's lensing galaxy is a large red elliptical, while 
the lensed source is blue, as expected for lens systems. 
The identified quad geometry is a typical le ns geome- 
try a nd may be called "cusp dominated" (e.g. iKochaneS 
l2QQ6f ). as it occurs when the source lies near to a cusp 
of the inner caustic. Although we also noted a strong 
disk component remaining after the initial Moffat profile 
subtraction, the robot is able to effectively model Danny 
as an SIS+external shear (a situation in which such a 
geometry would occur). 

Sandy has two strong arcs on either side of the lens 
galaxy, consistent with a double pattern produced by a 
source almost directly behind a SIS; this is well modeled 
by the robot. We note a fainter peak in surface brightness 
above the lens galaxy; with no obvious counter image, we 
believe this object is likely a small satellite galaxy in the 
lens plane. 

The two images comprising Kenickie's lensing event, 
and the center of the lensing galaxy, are not all perfectly 
aligned: this suggests that either the source does not lie 
quite on the optical axis of an ellipsoidal lens or that 



external shear is present. We note that the inner image 
appears to have more curvature than the outer; this could 
be brought about by unsubtracted lens-plane structure 
that a color image would rule out. 

Frenchie, another double, has been imaged in two fil- 
ters and has colors consistent with lensing. This is the 
only one of our systems backed by spectroscopic data. 
We find a spectroscopic redshift of Zd = 0.104. Also, a 
stellar velocity dispersion is available from SDSS, corre- 
sponding to (J* = 248 ±17 km s~^ after aperture correc- 
tion, in good agreement with the value estimated via the 
FP technique (226 ± 27 km s-^) 

Rizzo is able to be modeled as a double; however the 
morphologies of the two identifiable images are not as 
well-matched as in the previous 4 cases. In the lensing 
scenario, the source would lie only partly within the outer 
caustic of the lens, with the more extended outer image 
being only partially strongly-lensed. 

6.2. Lens environments 

We expect lenses, as massive galaxies, to be clustered. 
Consistent with this hypothesis, 3/4 new lenses are found 
at redshifts very similar to those of the main SLACS lens 
in the field. In this section we study the environment of 
the fields with more than one lens to identify possible 
large scale structures, using the environment measures 
of local and global density as defined in SLACS- VIII. 

SLACS- VIII found no significant difference between 
the environment of SLACS lenses and that of non- 
lensing, but otherwise identical, galaxies. With the ex- 
ception of Frenchie, the environments of the newly dis- 
covered lenses are fairly typical of the overall distri- 
bution of SLACS lens environments, with the field of 
Kenickie being slightly under-dense, while that of Sandy 
being somewhat over-dense. Frenchie, as previously men- 
tioned, lies in a very over-dense environment, the densest 
of ah SLACS fields. 

In fact, both Frenchie and Sandy appear associated 
with known clusters. Frenchie lies within SDSS-C41035, 
at redshift 0.106, which places both Frenchie {zd = 0.104) 
and its main lens {zd = 0.106) as members. The clus- 
ter in Sandy's field, MaxBCGJ217.49493+41. 10435, is 
at redshift 0.270, and the main lens at Zd = 0.285 lies 
within the cluster. We have only photometric redshift 
for Sandy {zd = 0.32 ± 0.01), which places it slightly be- 
yond the extent of the cluster; however given the possible 
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Ho^2m i/r«'=0.1 i/°f*=l.l 



if2=l I ifr<*'=0.05 ff°P*=2.7 



Fig. 3. — Lenses and lens models for new systems created by the HAGGLeS robot. We show the galaxy cutout (column 1), the subtraction 
residual (column 2), the masked robot input image (column 3), and the reconstructed images from the minimum source required to produce 
the configur ation (column 4). T he robot uses Moffat profile subtractions by default; here we have improved upon this by using a B-spline 
subtraction ([Bolton et al.|[2006r ). Cutouts are 6 arcseconds on a side. 
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Fig. 4. — Source redshift estimated from the FP method vs. the 
corresponding spectroscopic source redshifts for the SLACS lenses. 




Fig. 5. — Source redshift posterior probability distribution 
function estimated for the newly discovered lenses using the FP 
method. 

systematic errors in the photometric redshift we deem it 
likely that Sandy is also a member of the field's cluster. 

In conclusion, the results suggest that lens fields as- 
sociate with galaxy over-densities are the most likely 
to present additional strong lens phenomena. Only 
12/70 = 17% of ah SLACS lenses are associated with 
known clusters, but 50% of the newly discovered ones 
are (2/5 when including the likely candidate). With the 
present data it is hard to disentangle the contribution to 
the boost in strong lens surface density by the enhanced 
surface density of deflectors. 

6.3. Comparison to SLACS lenses 

We anticipated discovering lensing events similar to, 
but fainter than, the SLACS main lenses; we flnd this to 
be largely true. In Figure [H we plot apparent mag- 
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Fig. 6. — Plot of apparent i-band magnitude vs. lens redshift for 
the SLACS main lenses and new systems. We find that four of the 
five new systems are less luminous than the main lenses. 

nitude against redshift for the main lenses and new sys- 
tems. We flnd that four of the new systems are less lu- 
minous than the SLACS lenses at the same redshift, but 
that Frenchie is comparable. We show the distribution 
of stellar velocity dispersions in Figure [2 new systems 
are shaded. 

To interpret these histograms we need to take into ac- 
count the redshift dependence of the properties of the 
SLACS main sample. At the lower redshifts, the sample 
is dominated by the more abundant, slightly less mas- 
sive galaxies. At higher redshifts, the flux limit of the 
SDSS spectroscopic database leaves only the most mas- 
sive galaxies. Keeping this in mind, we also plot the 
distributions for the SLACS sample in the same red- 
shift range of all the new lenses save Frenchie (> 0.26 
is chosen because of a natural break in the redshift dis- 
tribution; the choice of threshold does not influence our 
conclusions). Indeed, the newly identifled lenses are less 
massive than the SLACS main lenses when both samples 
are restricted to Zd > 0.26. The average a* for the two 
samples are respectively 243 ± 35 km s~^ and 299 ± 20 
km 

6.4. Lens abundance 

In Table [HI we compare lensing rates (or abundances, 
in degree"^) for four HST imaging surveys, including 
this one. The quoted uncertainties on the inferred rates 
delimit the Bayesian 68% confldence interval, assuming 
Poisson statistics and a uniform prior PDF. As noted 
previously, classiflcation systems across surveys are not 
consistent, thus one must impose similar criteria when 
comparing the results of different surveys; naturally there 
will be differences in opinion. In this work we do not rely 
on spectroscopic data to classify a lens as "deflnite" and 
instead require that lens morphology, surface brightness, 
and color (if available) are clearly identiflable with a typ- 
ical lens geometry. Due to the limited data available, we 
apply these criteria rigorously. With the goal of applying 
similar criteria to all surveys considered, we take as defl- 
nite lenses: the literature-conflrmed lenses for MDS, the 
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TABLE 7 

Environments of new lenses and candidates 



Short name 


SDSS field 


local overdensity 


global overdensity 


Cluster 


ZCl 


Danny 


J0737+3216 


1.55±0.50 


1.22±0.28 


None 




Frenchie 


J1143-0144 


79±25.5 


3.7±0.50 


SDSS-C41035 


0.106 


Kenickie 


J1213+6708 


0.48±0.16 


0.45±0.16 


None 




Sandy 


J1430+4105 


1.83±0.62 


1.22±0.34 


MaxBCGJ217.49493+41. 10435 


0.270 


Rizzo 


J1103+5322 


0.98±0.3 


0.75±0.18 


None 





TABLE 8 
Lensing rates for selected surveys 



15 



Survey 


Area 


A^lenses 


Lensing rate 




(deg2) 




(degree"^) 


MDS^ 


0.17 


2 




AEGIS2 


0.19 


3 


16+15 


COSMOS^ 


1.64 


20 


12+1 


This Work 


0.13 


4 




with COSMOS prior 




(predicts 1.6) 





Confidence intervals are Bayesian 68%, assuming Poisson statistics 
and a uniform prior on the lensing rate, except in the final row 
where the COSMOS rate is taken as th e mean of an exponentia l 
prior. Relevant citations are as follows: MRatnatunga et al.l (119991 ). 
2 IMoustakas eFaP (l2007l ). ^ IFaure et all (|2008h 
"unambiguous" candidates for AEGIS, the "best" candi- 
dates for COSMOS, and the H = 3 lenses for this work. 

The largest survey o f "blank" sky undertaken to date is 
the COSMOS survey (jFaure et alboOST ). whose findings 
imply a measured lensing rate of around 12 lenses per 
square degree. Since our data are of comparable depth 
(and in the same single filter), we adopt this as our fidu- 
cial value. Given this lensing rate, we might expect to 
find 1.61 lenses, instead of our observed 4. Assuming a 
uniform prior on the lensing rate we find that the inferred 
lensing rate from our survey is 30tf lenses/de gree . 

The uniform prior is somewhat unrealistic - it does not 
down-weight the occasional high lens yields that could 
arise as statistical flukes from the long-tailed Poisson 
likelihood. A maximally conservative approach would be 
to take the COSMOS rate as the mean of an exponen- 
tial prior; in this case, we infer a lensing rate of 18^5^ 
lenses/degree^. We still find a significantly higher lens 
abundance than seen in the COSMOS survey: with the 
COSMOS prior the probability that the lensing rate in 
the SLACS fields is greater than 12 degree-^ is 88%. Re- 
laxing the assumption that our fields are similar in nature 
to those in COSMOS and reverting to the uniform prior, 
we find that there is only a ~ 2% chance that the lensing 
rate in the SLACS fields is less than the COSMOS rate 
of 12 lenses/degree^. 

Although no definitive conclusions can be drawn due 
to the small numbers and variety of surveys, the re- 
sults from this project do support the hypothesis of 
iFassnacht et al.l (|2QQ6l ) that looking near known lenses 
increases the efficiency of finding new lenses. 

7. CONCLUSIONS 

We have performed a highly efficient search for gravi- 
tational lenses based purely on imaging data. Our search 
area consisted of 44 HST/ ACS fields, each centered on a 
SLACS "definite" lens, therefore exploiting the expected 
clustering of gravitational lenses in each field. We com- 
pared the purity, completeness, and investigator time 
for four different search methods. These methods are 
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Fig. 7. — Histogram of central velocity dispersions for main lenses 
and new systems; new systems are shaded. Main lens velocity 
dispersion were aperture-corrected from SDSS photometry, while 
we used the FP relation to calculate the values for new lenses. 

comprised of: (1) use of the output from the HAGGLeS 
robot, (2) inspection of BG subtraction residual images, 
(3) inspection of BG cutout images, (4) full frame visual 
inspection of ACS fields. 
Our main conclusions are: 

• Taking into account efficiency as well as complete- 
ness and purity, we find that of the methods used, 
procedures 1 (using output from the HAGGLeS 
robot) and 2 (looking at subtraction residuals of 
bright galaxies) have the best performance. In situ- 
ations where the simple SIS+external shear model 
used by the robot may be insufficient - such as 
in clusters or to find naked cusp configurations - 
it may be most effective to use a combination of 
procedures 1 and 2, in order to inspect all objects 
modeled by the robot. However, looking at only 
the bright galaxies (BGs) in the fields did not de- 
crease the completeness of this particular survey, 
while doing so greatly improved efficiency. 

• We discovered four new strong lenses and one 
promising candidate in the course of our survey. 
We find that 3/4 of these new lens systems have 
lens redshifts similar to those of their main lenses; 
additionally, two of the new lenses are found in 
clusters of galaxies that also include their respec- 
tive main lenses. 

• We find that 3/4 new systems are less luminous 
and less massive than the SLACS lenses. Overall, 
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we are probing early-type galaxies at higher red- 
shifts and lower masses than the SLAGS survey. 
For these comparisons, we used the data from the 
improved robot models, available photometry, and 
the Fundamental Plane to estimate the central ve- 
locity dispersions and source redshifts for each of 
our new systems. 

• The lens abundance for this survey, 30lg^ 
lenses/degree^ (uniform prior), is markedly higher 
than the lensing rates for the three other HST sur- 
veys considered at comparable depths and resolu- 
tion. Despite the small numbers and variations in 
search methods, this result supports the idea that 
searching near known lenses increases the yield of 
a lens survey. 

The HAGGLeS project is currently using a combina- 
tion of procedures 1 and 2 to search for new lenses in 
the HST archive; through efforts such as this and oth- 
ers, we will be able to refine our lens search techniques 
for future surveys covering much larger areas of the sky. 
The use of efficient and repeatable lens search methods 
will further us towards the goal of having a large, ho- 
mogenous sample of strong gravitational lenses. Such a 
sample will enable us to calculate a lens-lens correlation 
function and constrain the statistical properties of halos 
containing lens galaxies. 
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ington. 
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Fig. 8. — Objects classified by the robot as if^^"^ > 2.5, including four main lenses (first and second rows). Cutouts are 6 arcseconds on 
a side. We find common false identifications to be rings left by the Moffat subtraction, unmasked disks, and nearby galaxies. 

APPENDIX 
SUPPLEMENTARY FIGURES AND TABLES 

In this Appendix we give a more complete illustration of the HAGGLeS robot's performance when ^iven the S LAGS 
main lenses. In Figure [8] we show all objects classified by the realistic robot as H^^^^ > 2.5: [Marshall et al.l (|2QQ8l ) 
found this to give a sample with ~ 20% completeness but ~ 100% purity. In the SLAGS fields, 10 objects were 
classified as H^^^^ > 2.5, including 4 main lenses. Gomparing to the human-classified results using the same input 
data (Procedure 2, Table [3|), this represents a completeness of ~ 11% and a purity of 40%. Some explanation for 
these differences are given in the main text. 

The classification of the SLAGS main lenses is listed in full in Table [9l The survey cutout images of these systems, 
sorted into bins in are shown in Figure[9](i^7^ > 2.5), Figure[l0](1.5 < H^p^ < 2.5), and Figure [l2](i^7^ < 1.5). 

For each system, we show the full cutout image as presented for inspection in Procedure 3, the lens galaxy-subtracted 
cutout image as presented for inspection in Procedure 2, and the lensed images and counter-images predicted by the 
HAGGLeS robot's best lens model. 
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Fig. 9. — Robot output for SLACS main lenses with if^^* > 2.5, including: galaxy cutouts (col. 1,4), residuals after Moffat subtraction 
(col. 2,5), and predicted images from the robot lens models (col. 3,6). In this figure and the two following, cutouts are 6 arcsec on a side. 
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Fig. 10. — Robot output for SLACS main lenses with 1.5 < Hr^ < 2.5, including galaxy cutouts, subtraction residuals, and predicted 
images from the robot lens models. 
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Fig. 11. — Continuation of Figure [TO] 
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Fig. 12. — Robot output for SLACS main lenses with Hr^ < 1.5, including galaxy cutouts, subtraction residuals, and predicted images 
from the robot lens models. 



